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The  total  momentum  of  a  thermodynamically  closed  system  is  unique,  as  is  the  total  energy. 
Nevertheless,  there  is  continuing  confusion  concerning  the  correct  form  of  the  momentum  and  the 
energy-momentum  tensor  for  an  electromagnetic  field  interacting  with  a  linear  dielectric  medium. 
Here  we  investigate  the  energy  and  momentum  in  a  closed  system  composed  of  a  propagating  elec¬ 
tromagnetic  field  and  a  negligibly  reflecting  dielectric.  The  Gordon  momentum  is  easily  identified 
as  the  total  momentum  by  the  fact  that  it  is,  by  virtue  of  being  invariant  in  time,  conserved.  We 
construct  continuity  equations  for  the  energy  and  the  Gordon  momentum  and  use  the  continu¬ 
ity  equations  to  construct  an  array  that  has  the  properties  of  a  traceless,  diagonally  symmetric 
energy-momentum  tensor.  Then  the  century-old  Abraham-Minkowski  momentum  controversy  can 
be  viewed  as  a  consequence  of  attempting  to  construct  an  energy-momentum  tensor  from  continuity 
equations  that  contain  densities  that  correspond  to  nonconserved  quantities. 


I.  INTRODUCTION 

The  energy-momentum  tensor  is  a  concise  way  to  rep¬ 
resent  the  conservation  properties  of  the  flow  field.  For 
most  types  of  simple  flows,  the  energy- momentum  tensor 
is  well-defined,  with  the  notable  exception  of  the  elec¬ 
tromagnetic  field  in  a  linear  dielectric  material.  The 
Abraham-Minkoswski  controversy  [Ir9]  for  the  momen¬ 
tum  of  electromagnetic  fields  in  a  dielectric  began  with 
the  derivation  of  the  energy-momentum  four-tensor  by 
Minkowski  0-  Noting  that  the  Minkowski  tensor  is  not 
diagonally  symmetric  and  therefore  contains  problems 
with  conservation  of  angular  momentum,  Abraham  0 
proposed  an  energy-momentum  tensor  that  was  symmet¬ 
ric,  but  at  the  expense  of  a  new  phenomenological  force. 
In  order  to  address  this  constraint  and  additional  issues, 
Einstein  and  Laub  13] ,  Nelson  and  others  proposed 
variants  of  the  energy-momentum  tensor. 

The  crux  of  the  Abraham-Minkowski  controversy  is 
whether  the  electromagnetic  momentum  density  in  a  di¬ 
electric  is  of  the  Minkowski  form 

Sm  =  —  (D  x  B)  (1.1) 

c 

or  the  Abraham  form 

gA  =  -(E  x  H).  (1-2) 

c 

Experimental  efforts  to  resolve  the  theoretical  impasse  in 
one  direction  or  the  other  have  not  been  definitive.  While 
some  experiments  favor  the  Abraham  formula,  other  ex¬ 
periments  support  Minkowski’s  version.  Brevick’s  0 
analysis  of  experiments  performed  by  Jones  and  Richards 
E3,  Ashkin  and  Dziedzic  [l4|,  and  others  showed  that 
the  allocation  of  momentum  between  the  field  and  ma¬ 
terial  was  the  determining  factor  in  whether  a  particular 
experimental  result  was  described  by  the  Abraham  or 
Minkowski  form  of  electromagnetic  momentum.  Follow¬ 
ing  Brevick  flij| ,  the  formula  for  the  field  momentum  has 
been  shown  repeatedly  to  be  essentially  arbitrary  such 
that  any  of  the  formulas  for  the  field  momentum  can  be 


combined  with  an  appropriate  momentum  for  the  mate¬ 
rial  to  produce  the  same  total  momentum  IH- 

In  1973,  Gordon  jlj  constructed  the  total  momentum 
from  a  microscopic  model  in  which  the  electromagnetic 
field  component  of  the  total  momentum  is  said  to  be  the 
Abraham  momentum  and  the  dielectric  is  treated  as  a  di¬ 
lute  collection  of  electric  dipoles  with  center-of-mass  mo¬ 
tion  in  the  direction  of  propagation  of  the  field.  Gordon 
[I]  discusses  the  empirical  and  experimental  validation  of 
the  total  momentum  density 

71 

gG  =  -(ExB)  (1.3) 

and  shows  that  the  density  gG,  integrated  over  a  volume 
containing  the  entire  field,  is  invariant  in  time.  Mikura 
flhj  subsequently  obtained  the  total  momentum  using  a 
more  general  material  model  that  included  polarizabil¬ 
ity,  magnetizability,  electrostriction,  and  other  material 
considerations.  Since  then,  the  total  momentum  has  been 
constructed  using  a  variety  of  field  and  material  momenta 
0- 

The  historical  development  of  the  electromagnetic 
energy-momentum  tensor  for  the  held  in  a  dielectric  has 
been  intertwined  with  efforts  to  derive  the  correct  elec¬ 
tromagnetic  momentum  density.  The  original  Minkowski 
and  Abraham  tensors,  like  the  Minkowski  and  Abraham 
momenta,  are  now  regarded  as  being  associated  with 
the  held,  alone,  requiring  a  material  component  to  com¬ 
plete  the  total  energy-momentum  tensor  [l5| .  Although 
a  number  of  composite  energy-momentum  tensors  have 
been  constructed  0  ,  the  tensor  forms  that  have  been 
constructed  from  held  and  material  constituents  fail  to 
satisfy  the  constraints  imposed  on  the  symmetry,  trace, 
or  four-divergence  of  a  total  energy-momentum  tensor. 

The  basic  tenant  of  classical  continuum  electrodynam¬ 
ics  is  that  the  electrodynamic  properties  of  a  linear  ma¬ 
terial  can  be  characterized  by  macroscopic  parameters. 
In  this  long  wavelength  limit,  the  material  is  modeled  as 
a  simple  linear  dielectric  with  a  linear  refractive  index 
n.  In  reflection,  the  electromagnetic  held  exerts  radia¬ 
tion  pressure  on  the  dielectric.  Neither  the  motion  of  an 
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unrestrained  dielectric  nor  the  forces  of  restraint  on  a  sta¬ 
tionary  dielectric  are  coupled  to  the  macroscopic  Maxwell 
equations  of  motion  for  the  electromagnetic  field.  For 
that  reason,  we  limit  our  consideration  to  the  case  of  neg¬ 
ligible  reflection  in  which  the  consequences  of  radiation 
pressure  vanish.  The  total  electromagnetic  momentum 
for  the  closed  system  consisting  of  the  macroscopic  elec¬ 
tromagnetic  field  and  a  negligibly  reflecting  linear  dielec¬ 
tric  is  unique.  Here,  we  identify  the  Gordon  momentum 
[1]  as  the  total  momentum  by  the  fact  that  it  is  invariant 
in  time  and  therefore  conserved  in  an  isolated  system. 
The  total  energy  is  shown  to  be  similarly  conserved.  In 
order  to  satisfy  the  constraints  imposed  on  an  energy- 
momentum  tensor,  we  perform  a  simple  change  of  vari¬ 
ables  and  derive  continuity  equations  for  the  total  energy 
and  the  total  momentum.  The  continuity  equations  are 
then  used  to  construct  an  array  that  has  the  properties 
of  a  traceless,  diagonally  symmetric  energy-momentum 
tensor,  but  in  a  coordinate  system  with  time-like  coordi¬ 
nate  ct/n. 

II.  ENERGY-MOMENTUM  TENSOR  OF 
NONINTERACTING  PARTICLES 


equations 


dpx 

-jrr  +  V  ■  p*u  =  0 
at 

(2.3a) 

¥+v ',a"=° 

(2.3b) 

¥+v  "‘u=a 

(2.3c) 

At  this  point,  we  adopt  a  four-dimensional  notation 
where  repeated  indices  are  summed.  We  take  Roman  in¬ 
dices  to  run  from  1  to  3  and  we  identify  the  coordinates 
xl  with  the  Cartesian  coordinates,  such  that  x1  =  x , 
x2  =  y ,  and  x3  =  z.  Greek  indices  run  from  0  to  3 
and  x°  is  identified  with  the  time-like  coordinate  ct.  The 
Minkowski  space-time  metric  is  diag(—  1, 1, 1, 1).  Finally, 
partial  differentiation  with  respect  to  the  indexed  coor¬ 
dinates  is  represented  by  da  =  d/(dxa). 

The  four  continuity  equations  m  and  m  can  be 
concisely  represented  by 

dpG 010  =  0,  (2.4) 


In  the  continuum  limit,  the  density  of  any  property  of 
identical  noninteracting  particles  can  be  treated  as  the 
number  density  multiplied  by  the  amount  of  the  property 
that  can  be  attributed  to  each  particle.  The  continuity 
equation  corresponding  to  a  specific  property,  such  as 
mass,  charge,  or  energy,  is  then  obtained  by  substitution 
of  the  specific  property  density  for  a  placeholder  number 
density.  For  an  infinitesimal  element  of  volume  in  an 
inviscid  sourceless  flow,  the  continuity  equation 

|f+V-pu  =  0  (2.1) 

is  derived  by  applying  the  divergence  theorem  to  a  Taylor 
series  expansion  of  the  property  density  field  p  and  the 
vector  velocity  field  u  =  ( ux ,  uy ,  uz)  of  the  flow  H3-  The 
continuity  equation  reflects  the  conservation  of  a  contin¬ 
uous  scalar  property  in  a  flow  in  terms  of  the  equality  of 
the  net  rate  of  flux  out  of  the  volume  and  the  time  rate 
of  change  of  the  property  density  p  inside  the  volume. 
Depending  on  the  context,  the  three- vector 


where 


i-  2 

czp 

cux  p 

cuyp 

cur  p- 

cpx 

uxpx 

uvpx 

uzpx 

cpv 

uxpy 

uypy 

uzpy 

cpz 

uxpz 

uypz 

uzpz. 

In  order  for  this  representation  to  be  unique,  we  must 
specify  the  relation  between  the  continuity  equations, 
that  is,  the  relation  between  elements  in  different  rows 
of  in  Eq.  (ELS  [3.  For  a  closed  system,  conserva¬ 
tion  of  angular  momentum  requires  Ga 13  to  be  symmetric 

0 

Ga' 3  =  Gfia.  (2.6) 

Therefore,  the  vector  property  p  must  represent  the  flux 
of  the  conserved  scalar  quantity  p  such  that  p  =  pu.  Con¬ 
versely,  the  flux  pu  of  a  scalar  property  of  the  particles 
is  a  conserved  vector  property  of  the  flow.  The  covariant 
form  of  the  continuity  equation  is  the  four-divergence, 
Eq.  (|2.4[).  of  the  energy-momentum  four-tensor 


g  =  pu  (2.2) 

is  known  as  the  momentum  density,  the  flux,  or  the  cur¬ 
rent  density  of  the  property.  Specifically,  g  corresponds 
to  the  linear  momentum  density  if  p  is  a  mass  density 
and  corresponds  to  the  charge  current  density  if  p  is  the 
electric  charge  density. 

Some  conserved  properties,  such  as  momentum,  are 
vectors.  For  a  flow,  the  density  of  a  conserved  vector 
property  can  be  represented  as  p  =  (px ,  pv ,  pz).  Applying 
the  scalar  formalism  to  the  three  orthogonal  components 
of  the  property  density  vector  yields  scalar  continuity 


"  c2p 

cpux 

cpuv 

cpuz  " 

cpux 

puxux 

P'UXUy 

puxuz 

cpuy 

puyux 

puyuy 

puyuz 

cpuz 

puzux 

puzuv 

puzuz_ 

The  continuity  equation  (12.41)  with  the  four-tensor  (12.71) 
is  valid  for  any  conserved  extensive  quantity  in  a  simple 
flow,  not  just  for  the  mass  of  a  fluid. 

The  energy-momentum  tensor  (ED  has  some  essential 
properties.  First,  the  four-divergence  of  each  row  vector 

dgGafj  =  0 


(2.8) 
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is  a  continuity  law  corresponding  to  the  conservation  of 
the  property  represented  by  the  property  density  in  the 
first  element  of  that  row.  Second,  the  tensor  is  diagonally 
symmetric 

Gal 3  =  GPa  (2.9) 

corresponding  to  the  absence  of  unbalanced  shear  forces 
and  conservation  of  angular  momentum  in  a  closed  sys¬ 
tem  0.  Third,  as  a  consequence  of  diagonal  symmetry, 

daGa P  =  0,  (2.10) 

the  four-divergence  of  each  column  vector  is  a  continuity 
law  corresponding  to  the  conservation  of  the  property 
represented  by  the  property  density  in  the  first  element 
of  that  column.  It  should  also  be  noted  that  the  simple 
energy-momentum  tensor  is  based  on  the  properties  of 
an  unimpeded  flow.  If  the  flow  is  redirected  by  impact 
with  a  macroscopic  object  then  one  is  obligated  to  include 
the  equations  of  motion  of  the  object  or  the  forces  of 
restraint. 


III.  THE  ABRAHAM  AND  MINKOWSKI 
ENERGY-MOMENTUM  TENSORS 


The  Abraham  and  Minkowski  energy-momentum  ten¬ 
sors  are  examples  of  a  number  of  different  tensors  that 
have  been  been  proposed  for  the  electromagnetic  field  in  a 
dielectric  0] .  The  Minkowski  tensor  can  be  constructed 
from  the  continuity  equations  for  energy  flux  and  mo¬ 
mentum  flux  by  the  same  procedure  that  was  used  to 
construct  the  array  (ESI).  The  Minkowski  tensor  is  not 
diagonally  symmetric  and  therefore  would  admit  viola¬ 
tions  of  angular  momentum  conservation  if  it  were  to 
be  considered  the  total  energy-momentum  tensor  of  a 
closed  system.  Instead,  the  Minkowski  tensor  is  consid¬ 
ered  to  be  a  representation  of  the  energy  and  momen¬ 
tum  of  the  field-only  component  of  a  system.  In  this 
section,  we  outline  the  construction  of  the  Minkoswski 
tensor  and  describe  how  the  procedure  is  modified  to  ob¬ 
tain  the  Abraham  tensor.  Neither  the  Minkowski  tensor, 
nor  the  Abraham  tensor,  satisfy  the  requirements  of  a 
total  energy-momentum  tensor.  The  Minkowski  tensor 
is  not  symmetric  while  the  Abraham  tensor  contains  a 
phenomenological  volume  force. 

The  macroscopic  Maxwell  equations  of  continuum  elec¬ 
trodynamics  are  the  basis  for  deriving  continuity  equa¬ 
tions  for  electric  and  magnetic  fields  in  a  dielectric.  For 
a  dielectric  with  no  free  charges  in  a  regime  of  negligible 
absorption  and  dispersion,  the  Maxwell  equations  may 
be  written  as 

<9B 

VxE  =  -9R  (3-la) 


V  x  B  = 


dn2  E 
d(ct) 


(3.1b) 


V  -B  =  0 


V-  ^E  =  0 


(3.1c) 

(3-ld) 


in  Heaviside-Lorentz  units.  The  electric  and  magnetic 
fields  can  be  defined  in  terms  of  the  vector  potential  A 
as 


dA 

d(ct) 

(3.2a) 

B  =  V  x  A 

(3.2b) 

for  transverse  fields  in  the  Coulomb  gauge. 

The  macroscopic  Maxwell  equations  are  the  axioms  of 
classical  continuum  electrodynamics.  Poynting’s  theo¬ 
rem, 


d 

d(ct) 


i  (n2E2  +  B2) 


+  V  •  (E  x  B)  =  0, 


(3.3) 


can  be  derived  by  multiplying  the  Faraday  law  (13.  lal) 
by  B  and  adding  it  to  the  Maxwell  -Ampere  law  (I3.1bl) 
multiplied  by  E.  Poynting’s  theorem  can  also  be  derived 
by  substituting  the  Maxwell  equations  into  the  temporal 
derivative  of  the  energy  density 


pe  =  (l/2)(n2E2  +  B2)  (3.4) 


using  a  vector  triple-product  identity.  Poynting’s  theo¬ 
rem  is  a  continuity  equation  for  the  Poynting  energy-flux 
vector 


SP  =  c(E  x  B)  =  (s^,4,4).  (3.5) 


The  theorem 


d 

d(ct) 


(?r2E  x  B) 


-B  x  (V  x  B)  +  B(V  •  B) 


-  n2E  x  (V  x  E)  +  E(V  •  n2E)  (3.6) 


is  derived,  in  a  manner  similar  to  Poynting’s  theorem, 
by  substituting  Maxwell’s  equations  into  the  temporal 
derivative  of  the  Minkowski  momentum  density,  g m  = 
(n2/c)(E  x  B).  The  right-hand  side  of  Eq.  (13.61)  can  be 
recast,  approximately,  as  the  negative  of  the  divergence 
of  the  Maxwell  stress  tensor  [0  with  components 


Wij  = 


-n2EiEj  —  BtBj  +  -  (n2E  -E  +  B-B)  <% 


(3.7) 

where  terms  involving  the  gradient  of  n2  have  been  ne¬ 
glected.  Then  the  temporal  derivative  of  the  Minkowski 
momentum  density,  Eq.  (EH,  can  be  expressed  using  the 
vector  divergence  operator,  V  •,  as 


—  (n2ExB)+V-W  =  0. 
o(ct) 


(3.8) 
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The  array 


'  Pe 

Sp/c 

Sp/c 

Sp/c 

c9m 

w11 

W12 

W13 

c9m 

W21 

w22 

w23 

-c9m 

W31 

w32 

w33_ 

(3.9) 


known  as  the  Minkowski  energy-momentum  tensor,  is 
constructed  from  continuity  equations  m  and  m- 
Using  the  summation  convention,  we  can  write  Poynt- 
ing’s  theorem  in  Eq.  (13.311  as 

dolM+djI^  =  0  (3.10) 

for  continuity  of  the  energy  flux  and  Eq.  (13.81)  as 

d0T^+djTZ  =  0  (3.11) 


for  continuity  of  the  momentum  flux.  Then,  each  row  of 
Tm  corresponds  to  a  four- divergence 

dpTl?  =  0.  (3.12) 

Because  the  array  T^f  in  Eq.  (13.91)  is  not  diagonally  sym¬ 
metric,  Abraham  proposed  the  energy-momentum  tensor 


'  Pe 

Sp/c 

sp/c 

s3p/c\ 

c9a 

W11 

W12 

w13 

c9a 

w21 

w22 

w 23 

-c9a 

w31 

w32 

w33_ 

(3.13) 


where  g.4  is  given  by  Eq.  m-  The  four  divergence  of 
this  tensor  is 


dpTf  =  -fp,  (3.14) 

where  /;3  is  the  Abraham  force.  The  Cartesian  compo¬ 
nents  of  the  Abraham  force 

i  =  WT)  (("2-1)ExB)  (3-15) 

are  obtained  by  substituting  the  Abraham  momentum 
density  into  the  continuity  equation  (13.81)  and  the  time- 
like  coordinate  component  is  /°  =  0. 

It  has  been  widely  reported  in  the  literature  that  nei¬ 
ther  Tm  nor  Ta  is  the  total  energy-momentum  tensor 
0-  Instead,  they  are  considered  to  be  two  of  many  ar¬ 
bitrary  forms  of  the  electromagnetic  part,  Tfld,  of  a  total 
energy  momentum  tensor 


T  =  Tfld  +  Tmati  (3.16) 

composed  of  energy-momentum  tensors  for  the  field  and 
material  subspaces,  however  those  subspaces  are  defined. 


IV.  TOTAL  MOMENTUM 

A  well-defined  quantity  for  the  physical  momentum  is 
not  derivable  solely  from  the  macroscopic  Maxwell  equa¬ 
tions  without  additional  conditions.  However,  the  total 


momentum  of  an  isolated  system  must  be  constant  in 
time.  In  this  section,  we  identify  the  unique  total  mo¬ 
mentum  using  this  constraint.  We  consider  the  case  of  a 
quasi-monochromatic  electromagnetic  field,  in  the  plane- 
wave  limit,  entering  a  linear  medium  from  vacuum  at 
normal  incidence.  The  medium  is  taken  to  be  a  sim¬ 
ple  linear  dielectric  in  the  regime  of  negligible  dispersion 
and  absorption.  A  dielectric  slab  is  partially  reflecting 
and  the  material  acquires  a  kinematic  momentum  that  is 
twice  the  momentum  of  the  reflected  field,  but  in  the  di¬ 
rection  of  the  incident  field.  In  the  absence  of  reflection, 
there  is  no  momentum  given  to  the  dielectric  slab  and  it 
remains  stationary.  We  adopt  this  case  of  a  stationary 
dielectric  in  which  reflections  can  be  neglected  by  assum¬ 
ing  that  an  antireflection  coating  has  been  applied  to  the 
dielectric  or  that  the  refractive  index  of  the  dielectric  is 
only  slightly  greater  than  unity. 

Gordon  [l|  used  a  microscopic  model  of  the  dielectric 
as  a  vapor  of  weakly  polarizable  atoms  and  derived  the 
material  momentum  as  the  continuum  average  of  the  me¬ 
chanical  momentum  of  the  atoms.  Assuming  a  rarefied 
vapor  of  atoms  in  order  for  reflections  to  be  negligible, 
Gordon  obtained  the  momentum 

Gg=  [  dv  gG=  f  dv-( ExB)  (4.1) 

JV  Jv  c 

by  adding  the  material  momentum  to  the  Abraham  mo¬ 
mentum.  In  continuum  electrodynamics,  the  electrody¬ 
namic  properties  of  a  material  are  characterized  only  by 
a  macroscopic  refractive  index.  Therefore,  the  micro¬ 
scopic  origin  of  the  Gordon  momentum  ED  is  of  no 
consequence  in  the  formalism  of  continuum  electrody¬ 
namics.  Following  Gordon,  we  will  demonstrate  that  the 
momentum  siD  is  invariant  in  time  and,  because  it  is 
conserved,  can  be  identified  as  the  total  momentum.  Be¬ 
cause  the  Gordon  momentum  depends  on  the  material 
only  through  the  refractive  index,  it  is  the  unique  total 
momentum  for  all  cases  in  which  the  medium  behaves, 
to  a  good  approximation,  as  a  linear  dielectric  with  re¬ 
fractive  index  n. 

Propagation  of  a  field  in  a  linear  medium  is  governed 
by  the  wave  equation, 


V2A 


n2  <92A 
~^~d¥  =  0, 


(4.2) 


written  in  terms  of  the  vector  potential  A,  B  =  V  x  A. 
For  quasi-monochromatic  plane  waves,  it  is  convenient 
to  write  the  vector  potential  in  terms  of  a  slowly  vary¬ 
ing  envelope  function,  A(z,t),  a  rapidly  varying  carrier, 
and  a  unit  vector,  e k,  in  the  direction  of  propagation  as 
A  =  A(z,t)e~dulot~kz'> e^.  For  a  plane-wave  entering  a 
dielectric  at  normal  incidence,  reflections  are  negligible  if 
Sn  =  n  —  1  is  small.  In  this  limit,  there  is  no  momentum 
given  to  the  bulk  material,  which  remains  stationary,  and 
we  can  apply  the  Fresnel  relation 


At  =  2 

Ai  n+1 


(4.3) 
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where  Ai  is  the  incident  amplitude  and  At  is  the  trans¬ 
mitted  amplitude.  Comparing  Eq.  USD  with  a  series 
expansion  of  1  /y/n  in  the  limit  of  small  Sn ,  we  find  that 
the  vector  potential  amplitude  inside  the  dielectric,  At, 
is  reduced  by  a  factor  of  y/n  from  the  incident  amplitude 
Ai. 

For  plane  cw  waves,  and  approximately  for  slowly  vary¬ 
ing  waves,  the  relation  between  the  amplitudes  of  the 
fields  simplifies  to 

|B|=n|E|  =  ™|A|.  (4.4) 

Then,  the  electromagnetic  energy  density,  pe  = 
(1/2)  (n2E2  +  B2),  can  be  written  as 

Pe  =  — |A|2.  (4.5) 

Applying  the  Fresnel  amplitudes  to  relate  the  fields  in¬ 
side  and  outside  the  medium,  At  =  Ai/ y/n ,  we  find  that 
the  energy  density  inside  the  material  is  a  factor  of  n 
larger  than  the  energy  density  of  the  same  field  in  the 
vacuum.  The  region  occupied  by  the  field  in  the  mate¬ 
rial  is  compressed  spatially  by  a  factor  of  n  due  to  the 
reduced  speed  of  light  within  the  medium,  such  that  the 
total  energy 

r  r  n2uj2 

U=  dvpe  =  /  dv — ^ —  I A|2  (4.6) 

Jv  Jv  c 

is  conserved.  Numerical  solutions  of  the  wave  equation 
for  a  held  entering  a  linear  material  through  a  gradient- 
index  anti-reflection  coating  indicate  that  the  held  in  the 
material  is  a  factor  of  y/n  smaller  and  a  factor  of  n  nar¬ 
rower  than  the  held  in  the  vacuum,  independent  of  the 
magnitude  of  ra,  as  long  as  reflections  are  suppressed  [20[  ■ 

Having  demonstrated  the  conservation  properties  of 
the  electromagnetic  energy,  we  demonstrate  the  conser¬ 
vation  properties  of  the  Gordon  electromagnetic  momen¬ 
tum  by  a  similar  procedure.  The  Gordon  momentum  [lj 
is  obtained  by  integrating  the  momentum  density 
over  all  three-dimensional  space.  Comparing  the  Gordon 
momentum,  expressed  in  terms  of  the  envelope  functions 

Gg  =  f  dv—  (E  x  B)  =  [  dv^\A\2ek,  (4.7) 
Jv  c  Jv  c 

with  the  total  energy,  Eq.  m,  we  see  that  conserva¬ 
tion  of  total  energy  implies  conservation  of  the  Gordon 
momentum.  The  total  momentum  of  a  closed  system  is 
unique  and  the  Gordon  form  of  total  momentum  is  con¬ 
served.  Therefore,  the  Gordon  momentum  can  be  iden¬ 
tified  as  the  total  momentum  of  the  thermodynamically 
closed  system.  By  uniqueness,  neither  the  Minkowski 
momentum  nor  the  Abraham  momentum  is  conserved. 

We  point  out  that,  in  the  macroscopic  limit  in  which 
the  dielectric  is  described  by  a  refractive  index  n,  the 
question  of  what  portion  of  the  energy  given  by  Eq.  ED 
resides  in  the  field  or  dielectric  is  improperly  posed.  Com¬ 
paring  Eqs.  gg  and  (1X71).  we  see  that  the  same  holds 
true  for  the  apportionment  of  the  momentum  into  field 
and  dielectric  components. 


V.  TOTAL  ENERGY— MOMENTUM  TENSOR 

In  the  previous  section,  we  identified  the  unique  to¬ 
tal  momentum  for  the  system  of  an  electromagnetic  field 
in  a  dielectric  for  the  case  of  negligible  reflections.  In 
this  section,  we  construct  the  corresponding  total  energy- 
momentum  tensor  from  continuity  equations  for  the  en¬ 
ergy  and  momentum. 

A  continuity  equation  is  a  differential  form  of  a  con¬ 
servation  law  applied  to  an  element  of  volume  in  a  con¬ 
tinuous  flow.  In  Sec.  II  we  showed  that  the  energy- 
momentum  tensor  for  dust  is  constructed  from  continu¬ 
ity  equations  in  which  the  differential  operators  act  on 
the  densities  of  conserved  quantities.  The  operand  of 
the  time  derivative  in  the  continuity  equation  (13.81)  is 
the  Minkowski  momentum  density.  Because  the  volume 
integral  of  the  Minkowski  momentum  density  is  not  a 
conserved  vector  quantity,  we  do  not  consider  Eq.  m 
to  be  a  suitable  continuity  equation  with  which  to  con¬ 
struct  an  energy-momentum  tensor.  Instead,  we  write 
Eq.  (13.81)  as 

TL  d 

--cgG+V-W  =  0,  (5.1) 

c  at 

in  terms  of  the  Gordon  momentum  density  gG  = 
(n/c)(E  x  B)  whose  volume  integral  is  a  conserved  vec¬ 
tor  quantity.  Equation  (E3D  provides  three  continuity 
equations  for  our  energy-momentum  tensor.  The  addi¬ 
tional  continuity  equation  is  obtained  by  writing  Poynt- 
ing’s  theorem  in  (13.31)  as 

+V  •  [n(E  x  B)]  =  —  -n(E  x  B)  (5.2) 
c  at  n 

using  the  densities  of  conserved  quantities.  We  are  con¬ 
sidering  the  case  of  a  closed  system  in  which  there  are 
no  reflections.  As  this  is  assumed  to  be  accomplished  by 
a  gradient-index  antireflection  coating,  we  can  drop  the 
term  containing  Vn  and  write  Poynting’s  theorem  in  Eq. 
(15.21)  as 

^|+v-[n(ExB)]=°.  (5.3) 

In  Maxwell’s  equations,  expressed  in  terms  of  3- 
vectors,  time  is  not  a  coordinate.  We  make  a  simple 
change  of  time  variable  to  t  =  t/n  and  write  the  conti¬ 
nuity  theorems  ED  and  ED  as 


1  d 

~vrc  gc  +V-W  =  0 

C  OT 

(5.4a) 

--p  +  V  ■  cgG  =  0, 
c  or 

(5.4b) 

where  the  Gordon  momentum  density  is  given  in  Eq. 
ED-  However  when  writing  Maxwell’s  equations  as  ten¬ 
sor  equations,  time  is  one  of  the  four  space-time  coordi¬ 
nates  and  we  define  the  time-like  coordinate 

x°  =  ct  =  —.  (5-5) 
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Then  the  four  scalar  continuity  equations,  Eqs.  (I5.4al) 
and  (1 5 . 4  b  |l .  can  be  written  concisely  as  a  single  equation, 
as  in  Section  2.  We  define  an  operator 

d*  =  (^o,dx,dy,dz)  ^5'6) 

and  an  array 


j-ia/3  _ 

i 

O'-1** 

c9g 

W1 2 3 4 5 6 7 8 9 10 11 

c9g 

W12 

c9g~ 

W13 

c9g 

W21 

w 22 

W23 

-c9g 

W31 

w32 

w33 

such  that 

dpT01^  =  0.  (5.8) 

The  array  (15.71)  has  a  number  of  notable  properties. 
The  array  is  diagonally  symmetric 

Tal 3  =  TPa  (5.9) 

and  has  a  vanishing  trace 

T“  =  0.  (5.10) 

The  operator  (15.61)  applied  to  the  rows  of  (15.71)  generates 
a  continuity  law  for  a  demonstrably  conserved  electro¬ 
magnetic  energy  or  momentum  property.  The  operator 
(15.61)  applied  to  the  columns 

BaTaf)  =  0  (5.11) 


VI.  CONCLUDING  REMARKS 

For  many  years,  the  Abraham-Minkowski  contro¬ 
versy  has  been  resolved  by  postulating  a  total  energy- 
momentum  tensor  that  is  comprised  of  separate  field  and 
matter  tensors,  an  approach  that  involves  assumptions 
about  the  behavior  of  matter  in  the  presence  of  an  elec¬ 
tromagnetic  field.  In  continuum  electrodynamics  the  in¬ 
teraction  of  the  field  and  matter  is  described  in  terms 
of  a  single  macroscopic  parameter,  the  refractive  index 
n.  We  showed  that  the  total  momentum  of  a  thermody¬ 
namically  closed  system  consisting  of  a  quasimonochro- 
matic  field  and  negligibly  reflecting  linear  dielectric  in 
the  continuum  limit  is  the  Gordon  momentum.  We  de¬ 
rived  continuity  equations  from  the  Maxwell  equations 
and  used  a  time  variable  transformation  to  write  the 
continuity  equations  in  terms  of  densities  of  conserved 
energy  and  Gordon  momentum  quantities.  When  writ¬ 
ten  in  four-dimensional  tensor  form  with  time-like  coor¬ 
dinate  ct/n ,  the  continuity  equations  are  obtained  from 
the  four-divergence  of  a  traceless,  diagonally  symmetric 
energy-momentum  tensor. 

In  summary,  we  constructed  the  energy-momentum 
tensor  in  Eq.  EZJ  from  continuity  equations  that  were 
derived  from  the  macroscopic  Maxwell’s  equations.  It  is 
interesting  to  note  that  a  time  coordinate  transformation 
was  required  in  order  to  write  the  continuity  equations  as 
the  four-divergence  of  the  symmetric  energy-momentum 
tensor. 


generates  the  same  continuity  equations  as  a  consequence 
of  symmetry.  These  are  the  properties  that  we  associate 
with  an  energy-momentum  tensor. 
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